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We consider the liquid-vapor type phase transition for fluids conflned within spatially periodic 
external fields. For a fiuid in d = 3 dimensions, the periodic field induces an additional phase, 
characterized by large density modulations along the field direction. At the triple point, all three 
phases (modulated, vapor, and liquid) coexist. At temperatures slightly above the triple point and 
for low (high) values of the chemical potential, two-phase coexistence between the modulated phase 
and the vapor (liquid) is observed. We study this phenomenon using computer simulations and 
mean-field theory for the Ising model. The theory shows that, in order for the modulated phase to 
arise, the field wavelength must exceed a threshold value. We also find an extremely low tension of 
the interface between the modulated phase and the vapor/liquid phases. The tension is of the order 
10"* feflT per squared lattice spacing, where fcs is the Boltzmann constant, and T the temperature. 
In order to detect such low tensions, a new simulation method is proposed. We also consider the case 
of d = 2 dimensions. The modulated phase then does not survive, leading to a radically different 
phase diagram. 

PACS numbers: 64.70.F-,64.75.-g,75.40.Mg 



I. INTRODUCTION 

Liquid-vapor type phase transitions in fluids are pro- 
foundly affected by confinement (for a recent review see 
Ref. [l|) . Typical effects are the depression of critical tem- 
peratures [J , changes in universality Q , or entirely new 
phenomena altogether 0, [H . The confinement of a fluid 
between two parallel surfaces is arguably the most sim- 
ple example one could envision Q . Already for this case 
the corresponding phase behavior is extremely rich, es- 
pecially if the surfaces have different interactions with 
the fluid Q . With the advance of microcontact print- 
ing 0, vapor deposition and grafting methods Q, as 
well as photolithography the possibilities of tuning 
the surface-fluid interaction are essentially endless. In 
addition to surfaces, confinement in fluids may also be 
induced via external fields (for example, optical tweezers 
can be used to realize one-dimensional diffusion channels 
for colloidal particles in suspension [I^l)- Hence, well- 
characterized geometries of ever increasing complexity 
can be generated, and the phase behavior of fluids con- 
fined within these is expected to become correspondingly 
richer. 

With these developments in mind, this paper considers 
the fate of the liquid-vapor transition in a fluid confined 
within a static external field having periodic spatial os- 
cillations in one direction. In c? = 2 dimensions, such a 
field might be realized using a stripe-patterned surface 
pdl while in d = 3 dimensions, laser [l^ or electric 
fields p^l - [la | could possibly be used. The case d = 3 
was first considered theoretically in Ref. [l3 for a colloid- 
polymer mixture. The main finding was a new kind of 
phase transition, referred to as laser-induced condensa- 



tion (Lie), which takes place provided the field wave- 
length is large enough. In the presence of the periodic 
field, one then observes a new third phase (in addition to 
the vapor and liquid phases) characterized by (i) an aver- 
age density between that of the vapor and liquid phase, 
and (ii) featuring large density modulations along the 
field direction (because of the latter modulations we re- 
fer to this phase as the "zebra" phase in what follows). 
The presence of the zebra phase dramatically alters the 
liquid- vapor phase diagram: the critical point of the bulk 
transition is replaced by two new critical points and a 
triple point. At temperatures between the triple and 
critical points, vapor-zebra and liquid-zebra two-phase 
coexistence is observed (at low and high values of the 
chemical potential, respectively). 

In a subsequent publication [isj the nature of the crit- 
ical points was elucidated, and also the tensions j^z and 
7iz of, respectively, the vapor-zebra and liquid-zebra in- 
terfaces were calculated. The main observations were a 
critical behavior corresponding to effectively d = 2 di- 
mensions (i.e. one below the system dimension), and ex- 
tremely low interfacial tensions. The latter were found, 
using density functional theory, to be at most j^z ^ 7iz 
10~^ ksT per projected particle area (with ks the Boltz- 
mann constant, and T the temperature). The accompa- 
nying simulations confirm that 7vz and jiz must be ex- 
tremely low, but no numerical values could be obtained 
(from the simulation data of Ref. [l^ interface tensions 
of exactly zero cannot be completely ruled out either). 

In this paper, we revisit LIC using computer simula- 
tions and mean-field theory for the Ising model. Com- 
pared to a colloid-polymer mixture, computer simula- 
tions of the Ising model allows for much faster equili- 
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bration, such that larger system sizes can be reached. In 
addition, the underlying spin reversal symmetry of the 
Ising model makes the finite-size scaling analysis much 
more straight-forward. Of course, since the universality 
class of fluids is the Ising one, generic trends observed 
in the latter directly apply to fluids as well. Wc first 
consider LIC in d = 3 dimensions. The corresponding 
phase diagram is calculated using both simulation and 
theory. In particular, we demonstrate how the phase di- 
agram depends on the field wavelength and amplitude. 
Our next aim is to measure the interface tensions 7vz 
and 7iz : it is important to confirm the density functional 
prediction that 7vz and 7iz are extremely low but finite. 
As it turns out, accurate measurements of such extremely 
low tensions are beyond the scope of "standard" methods 
[l^ll^, and so an alternative route is proposed. Finally, 
we consider LIC in c? = 2 dimensions. Since the critical 
behavior was shown to be of dimension d — 1 [isj , we 
expect radical departures from the previously considered 
case d = 3. Indeed, in c? = 2 dimensions, the two critical 
points do not survive, and an altogether different phase 
diagram is obtained. 



II. MODEL AND SIMULATION METHOD 

We consider the Ising model on rectangular L x L x D 
{d = 3) and LxD {d ^ 2) lattices with periodic boundary 
conditions in all directions. The system is exposed to a 
periodic external field ypcr(z), with the z-axes parallel 
to edge D of the lattice. To each lattice site i, a spin 
variable s.^ = ±1 is attached. The energy of the system 
is given by 

E = -JY, SiSj +H^s, + Y^ s,Vp^,{zi), (1) 

where the first sum is over nearest neighbors, and the 
remaining sums over sites. The first term is the usual 
Ising pair interaction with coupling constant J (we con- 
sider ferromagnetic interactions J > only) . The second 
term is the interaction of the spins with a homogeneous 
external magnetic field of strength H. The last term rep- 
resents the interaction with the periodic field, where Zi is 
the z-coordinate of spin i. For the periodic external field 
we use a block wave of alternating sign 



Vpcriz) 



-h 0<z< A/2, 
+h X/2 < z< A, 



(2) 



with h the field strength, and A the wavelength. Due to 
the discretization of the lattice we must choose A = 2ani, 
with a the lattice constant, and ni an integer. The use 
of periodic boundary conditions implies that the lattice 
edge D = Xn2, with n2 also an integer. In what follows, 
the lattice constant is the unit of length a = I. In addi- 
tion, a factor of l/fc^T is assumed to have been absorbed 
into the coupling constants J, H and h such that these 
quantities arc dimcnsionlcss. 



Monte Carlo simulations and mean-field theory are 
used to study the phase behavior of the above model. 
The key output of the simulations is the distribution 
P(m), defined as the probability of observing the sys- 
tem in a state with magnetization m = ii/N)'^^Si, 
with N = DU^^^ the total number of lattice sites. We 
emphasize that P{m) depends on all the model param- 
eters introduced above, including the system size. To 
obtain P{m) we use single spin-fiip dynamics [2lj com- 
bined with successive umbrella sampling j2^ : the latter 
scheme ensures that P{m) is obtained over the entire 
range — 1 < to < 1, including those regions where P(m) 
is very small. We also use histogram reweighting [23| 
to extrapolate data obtained for one set of values of the 
coupling constants to different (nearby) values. 



III. RESULTS IN rf = 3 DIMENSIONS 

In this section we present results for the case d = 3. 
We begin in Sec. IIII Al with simulation results obtained 
for an external potential with strength h = 0.075 and 
wavelength A = 10. Following this, in Sec. IIII Bl we 
present mean-field theory results for how the phase dia- 
gram varies as the parameters h and A are varied. 



A. Laser-induced condensation: phase diagram 

To understand LIC in the Ising model it is best to con- 
sider the free energy F(m) as function of the magnetiza- 
tion m. The latter is related to the magnetization prob- 
ability distribution, F{m) = —kBTlnP{m), to which we 
have direct access in our simulations. In Fig. [Ija), we 
show F{m) for a high value of the coupling constant J 
and H — 0. The salient features are two global min- 
ima, at low and high values of m, reflecting a coexistence 
between two phases (I and II). We also observe a local 
minimum at m = 0, corresponding to a phase III, but it 
is meta-stable. In (b), we plot F{m) for a lower value of 
J and H = 0. We now observe three minima at equal 
height corresponding to a triple point, where all three 
phases coexist. Next, in (c), we show F{m) for an even 
lower value of J, but still using H = 0. There is now only 
one global minimum at m = 0. However, by applying an 
appropriate homogeneous field _ff = Ai, a coexistence 
between phase I and III is obtained (d). The value of Ai 
follows from the slope of the "common tangent" line li. 
Similarly, by applying a homogeneous field H = A2 (de- 
termined from the slope of line I2), coexistence between 
phase II and III is obtained (e). Finally, at some critical 
value J — Jcr,i(^cr,2), thc I-III (II-III) coexistence line 
terminates, below which there is only one phase. 

Fig. [T] is the analogue of LIC [Ol in the Ising model, 
with phase I playing the role of vapor (v), phase II of 
the liquid (1), and phase III of the zebra (z) phase. Due 
to spin reversal symmetry, it holds that If = at the 
triple point J = Jtr- Below the triple point, symmetry 
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FIG. 1: Lie in the d — 3 Ising model; plotted in each of the graphs is the free energy F{m) in units of ksT (vertical axes) 
versus the magnetization m (horizontal axes). The free energy curves present actual simulation data obtained for system size 
L = 10 and D = X. (a) Free energy for J > Jti and H — 0, i.e. above the triple point. A coexistence between two phases, I 
and II, is observed, (b) Free energy for J — Jtr and H = 0, i.e. exactly at the triple point; three-phase coexistence is observed, 
(c) Free energy measured between the triple and critical points, < J < Jtr, but still using H = 0. We now observe two 
"common tangent" lines: Zi and h- By choosing H = Ai, where Ai is determined from the slope of h, coexistence between 
phases I and III can be induced (d). Similarly, from the slope of h, we obtain H = A2, at which phases II and III coexist (e). 



implies that Ai = — A2 = A and Jcr,i = Jcr,2 = Ja- 
The resulting phase diagram is a symmetric pitchfork 
(Fig.[2l^a)). The crucial difference with fluids (which typ- 
ically lack spin reversal symmetry) is that the phase di- 
agram is asymmetric in that case: Jcr,i 7^ Jcr,2 and the 
fields (chemical potentials) A^ are not trivially related to 
each other [isj . 

We emphasize that the free energy curves in Fig. [1] 
are obtained in simulations using D = X. If one in- 
stead uses D = Xn2 with integer ri2 > 1, one finds that 
F(m) develops additional minima, as discussed in detail 
in Ref. [l^ These additional minima refiect meta-stable 
coexistence states and should not be confused with new 
phases. Hence, also when D > X, the generic mechanism 
of Lie as shown in Fig. [T] still applies. 



B. Stability of the zebra phase: 
mean-field calculations 



In order to develop a qualitative understanding of how 
the Lie phase diagram of Fig. ^a) depends on the ex- 
ternal field wavelength A and amplitude h, we use the 
following simple mean-field (Bragg- Williams) approxima- 



tion [23, [23 for the free energy F of the system with 
Hamiltonian E, given in Eq. ([T|): 



+kBT 




(3) 



where rrii = (si) is the average magnetization at lat- 
tice site i. For a given external potential Vpcriz), the 
average magnetization profile corresponds to the set 
{mi, TO2, • ■ ■ } which minimize the free energy ([3]); i.e. are 
the solution to the set of equations dF/drrii = 0. This 
yields the following set of i simultaneous equations: 



In 



1 



1 — nij 



0, (4) 



where denotes the sum over the 6 (in = 3) nearest 
neighbor lattice sites of site i. Because the external po- 
tential Vpor(^) in Eq. only varies in the z-direction. 
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FIG. 2: (a) Phase diagram of LIC for the d = 3 Ising 
model. The phase diagram is a symmetric "pitchfork" fea- 
turing one triple point, and two critical points (indicated by 
dots). The lines correspond to first-order phase transitions 
(between phases I-II, I-III, and II-III). (b) LIC in d = 2 di- 
mensions. In this case, there is no phase III. The phase dia- 
gram features only a single line of first-order phase transitions 
terminating in a critical point. 
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FIG. 4: LIC phase diagrams for d — 3 obtained using mean- 
field theory (for clarity only the transition lines between the 
zebra III phase and the I and II phases are shown). In (a) 
we show results for h = 0.075 and various values of the field 
wavelength A. In (b) we show results for A = 10 and various 
values of the field amplitude h. 
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FIG. 3: A sequence of magnetization profiles for varying J, 
going from J = 0.14 (top) to J = 0.24 (bottom), in increments 
of 0.01 obtained from the mean-field theory. For low values 
of J, one observes phase III, which is characterized by large 
modulations in the magnetization profile. For high values of 
J, phase II is observed, which is characterized by an overall 
low value of the magnetization and only small modulations. 
Note the discontinuous change in the profiles as the system 
crosses the II-III phase boundary (here at J « 0.18 — 0.19). 
The remaining parameters used in the figure are H — 0.03, 
h = 0.075 and A = 10. 



we have magnetization profiles that only vary in this one 
direction and so solving Eqs. ^ is straightforward. We 
do so using a simple (Picard) iterative numerical scheme. 
In Fig. [3] we show some example magnetization profiles 
calculated for various different values of J as one crosses 
the transition line from phase II to the zebra phase III. 
We see a discontinuous change in the average magnetiza- 
tion in the system as one crosses the phase transition. 

In Fig. m^a) we show phase diagrams for various val- 
ues of the field wavelength A and fixed field amplitude 
h = 0.075 (the upper and lower horizontal lines corre- 
spond to zkh, respectively). Note that for clarity we only 
display the I-III and II-III coexistence lines and do not 
display the I-II liquid- vapor coexistence line. For A = 10, 
the zebra critical points arc marked A and i?, while point 
C indicates the triple point. In the limit X ^ oo, the 
critical points A and B shift toward ( Jcr,buik, re- 
spectively, where Jcr,buik = 1/6 in the mean- field theory. 
As A — >■ oo, essentially two infinite systems are obtained: 
one inside a positive (homogeneous) external field h, and 
one inside a negative field —h. The value of H at the 
respective critical point simply has to "cancel" this field. 
In the opposite limit A — )■ Amin = 2, we observe the loss 
of the zebra phase. In order for the zebra phase to sur- 
vive, A/2 must exceed the bulk correlation length which 
is the quantity that determines the distance over which 
the average density changes from one value to another. 
Indeed, for A = 4 and smaller, the critical points A and 
B can no longer be identified, and only point C survives 
(which then no longer is a triple point, but a critical 
point, marking the end of the I-II coexistence region). 
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When A = 4 and h = 0.075 the mean-field critical point 
is at J « 0.177 and when A = 2 it is at J « 0.169. Recall 
that the mean-field bulk critical point (i.e. for h ~ Q) is 
at Jcr.buik = 1/6 « 0.167. 

In Fig. m^b) we show phase diagrams for fixed A = 10 
(chosen above the threshold such that the zebra phase 
survives) and various values of the field amplitude h. In 
the limit /i — 0, we observe that the points A, B and 
C all approach ( Jcr.buik = 1/6, H = 0), the critical point 
of the bulk system. When h is very small, it is difficult 
to locate numerically the transition points. However, a 
threshold value of h below which the zebra phase vanishes 
appears to be absent in this case (in contrast to the case 
as A is decreased) . The effect of increasing h is that the I- 
III and II-III transition lines open-up, with the transition 
points A, B and C shifting toward larger values of J. 
Note that the value of H at the critical point is always 
less in magnitude than the value of h. When h = 0.02, 
then the critical value H^r ~ 0.005; when h ~ 0.2, then 
ifcr ~ 0.163 and when /i = 2, then H„ « 1.957. We see 
from these that as h becomes large, then _ffcr ^ h. 



C. Finite-size scaling analysis 

We now continue with our simulation analysis using 
A = 10 and h = 0.075. Finite-size scaling is used to 
locate the triple and critical points. We measure P{m) 
for various values of L, keeping D = X fixed. We thus 
assume that correlations in the z-direction are "cut-off" 
by the periodic field, and so we do not need to scale 
in this direction (we return to this point shortly). The 
distribution P{m) is always measured at H = and sym- 
metrized by hand afterward such that P{m) = P{—m), 
thereby imposing the spin reversal symmetry of the Ising 
model; subsequent histogram reweighting (in J and H) 
is performed using the symmetrized distribution. 

To determine Jtr we use H = and assume that, pre- 
cisely at the triple point, the magnetization distribution 
P(m) is a superposition of three (non-overlapping) Gaus- 
sian peaks, centered around m = —mo, m = 0, and 
771 = +mo, respectively. For such a triple-peaked dis- 
tribution one may easily show that the Binder cumulant 
Ui = 2/3, with the cumulant defined as 

r + l 

U4={mY/{m^), {m'') = J m''P{m)dm, (5) 

and where it is assumed that P{m) is normalized. In 
Fig.in^a), for the case when h = 0.075 and A = 10 we plot 
versus J for various system sizes L. The curves strik- 
ingly intersect at the expected "height" of a triple point, 
from which we conclude that Jtr ~ 0.2811. Note that 
this exceeds Jcr,buik ~ 0.2217 of the critical point in the 
bulk {h = 0) Ising model [2^ , consistent with the general 
observation that confinement lowers transition tempera- 
tures. In Fig. O^b), we plot the free energy F{m) at the 
triple point for various system sizes. The curves clearly 
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FIG. 5: Finite-size scaling analysis to locate the LIC triple 
point for the d = 3 Ising model, (a) Binder cumulant U4, 
versus the coupling constant J for various system sizes L. The 
curves strikingly intersect at f/4 — 2/3 (horizontal line) of a 
triple point; the value of J at the intersection yields Jtr. (b) 
The scaled free energy F{m)/L precisely at the triple point, 
J — Jtr, for various system sizes. Clearly visible are the three 
minima, corresponding to the coexisting phases, and a free 
energy barrier that increases linearly with L. 



show the three minima of the coexisting phases. Note 
that the minima are shifted to zero, and that the vertical 
scale is divided by L. In this representation, the barrier 
AF/L (vertical arrow) is approximately constant. Hence, 
at the triple point, we observe a free energy barrier that 
increases linearly with the system size AF ex L. This im- 
plies that the general shape of F{m), i.e. featuring three 
minima, persists in the thermodynamic limit L — > 00, 
and thus reflects a genuine triple point (see also Ref. [23 
where these ideas were first applied to first-order phase 
transitions). 

For values of J between the triple and critical points, 
coexistence with the zebra phase (phase III) is observed 
at appropriate values H = ±A of the external magnetic 
field. To locate Jcr we perform the same cumulant anal- 
ysis as in our previous work [l^. For a given value of 
J; (m) and are measured as function of H (due to 
symmetry, only H > needs to be considered). One 
then uses these data to construct a graph of U4 versus 
(?7i) (which thus is parametrized by H). The resulting 
curve reveals a maximurn, corresponding to _ff = A, en- 
veloped by two minima [2^. The average value of the 
cumulant at the minima equals Q4. In Fig.[6l^a), we plot 
Q4 versus J for different L; from the intersection point 
we conclude that Jcr ~ 0.2531 (the corresponding critical 
field Acr ~ 0.017). The difference M in the magnetiza- 
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FIG. 6: Finite-size scaling analysis to locate the LIC critical 
point for the d — 3 Ising model, (a) Cumulant Q4 versus the 
coupling constant J for various system sizes L; the intersec- 
tion yields J = Jcr. (b) Scaling plot of the order parameter 
M where t = J/ Jcr — 1; by using Jcr f« 0.2531 and d = 2 Ising 
values for the critical exponents /3, u the data for different L 
collapse onto a single master curve. 



tions (m) at the minima yields the order parameter. The 
latter is analyzed in the finite-size scaling plot [2l|, [2^ 
of Fig. ^h). The key result is that, by using the d = 2 
Ising values for the critical exponents {/3 = 1/8, ly = 1) 
the data for different L collapse. The critical points of 
LIC in d = 3 dimensions thus remain in the universality 
class of the Ising model, but of reduced dimensionality 
d — 1. We believe that colloid-polymer mixtures should 
ultimately yield the same exponents, but that their com- 
plexity still prevents efficient simulations of large enough 
systems to exphcitly see this [T8| . 

D. Measurement of ultra-low surface tension 

We now consider the regime between the triple and 
critical points using H = ±A. The free energy then 
schematically resembles Fig. [ljd,e), corresponding to T 
III and II-III phase coexistence, respectively. Hence, 
there will be interfaces present, and our aim is to measure 
the corresponding interface tension 7 (due to symmetry 
it holds that "fijii = "fiijii = 7: of course). Following 
density functional calculations [l8|, 7 is expected to be 
extremely low. In principle, for liquid-vapor transitions, 
the corresponding interface tension 71V can be accurately 
determined from the free energy using an idea of Binder 
[sot ■ I'^ '^^^ previous work [l^ , we discussed how this ap- 
proach may be generalized to LIC, but it was clear that 
present computer power is not sufficient to reach the sys- 
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FIG. 7: Instantaneous magnetization profiles m{z) obtained 
for J — 0.275, L = 14, and fixed overall magnetization 
?n = —0.426. In (a) we observe a coexistence between two 
domains, while in (b) a coexistence between four domains is 
seen. The vertical lines indicate the approximate locations 
of the interfaces. By counting how often the arrangements 
(a) and (b) occur during a long simulation run, the interface 
tension can be determined, see Eq. ((6)1. 



tem sizes required for this method to work. Hence, we 
propose a different method. 

The key observation is that the periodic field Vpcr{z) 
suppresses interface fluctuations (capillary waves) in the 
z-direction: even though 7 is very low, the I-III and II-III 
interfaces are sharp. This is in contrast to conventional 
liquid-vapor interfaces which, at low interface tension, 
are extremely broad [3l|, HH . The fact that the interfaces 
remain sharp is the property we exploit to extract 7. To 
this end, we consider a simulation box with edge D ~ 4A. 
In Fig. [71 we show instantaneous magnetization profiles 
m{z) obtained for two equilibrated samples at fixed over- 
all magnetization m. The value of m must be chosen such 
that half the system is occupied by phase I, and the re- 
mainder by phase HI, which can be obtained from the lo- 
cal maximum in the free energy (point "A" in Fig.[T](d)). 

Since the interfaces are essentially flat, one can eas- 
ily identify where the phases are located. In Fig. ^a) , 
we see one large domain of phase I (characterized by a 
low overall magnetization) coexisting with one large do- 
main of phase HI (characterized by large modulations 
in the magnetization). Hence, two I-III interfaces are 
present (recall that periodic boundaries are used). In 
Fig. [7{b), we again observe I-III phase coexistence, but 
this time the phases are arranged such that four I-III in- 
terfaces are present. In equilibrium, arrangement (a) is 
preferred since it has the smallest interface area: 2L^ ver- 
sus 4L^, with L the lateral box size. However, for finite 
L, arrangement (b) is also frequently observed, since 7 is 
small. In fact, from the ratio of counts R, the interface 
tension can be determined 

In i? = + A5, i? EE na/ub, (6) 

where denotes the number of times arrangement i = 



L 20 




FIG. 8: (a) Variation of In R with system size L for J — 0.275; 
symbols are simulation data, the curve is a fit to Eq. ((6]). 
Note that the data are consistent with lim_L_>oln-R = In 2 of 
the entropy difference, (b) The interfacial tension 7, in units 
of fcflT per squared lattice spacing, as a function of J. In 
agreement with theoretical expectations [l^, 7 is extremely 
low and decreases as J — > Jcr- 



(a, b) was seen during a very long simulation run (note 
that this simulation must be performed at fixed m cho- 
sen to yield equal volumes of both phases). The "off- 
set" AS reflects the combinatorial and translational en- 
tropy difFcrcncc between the arrangements. The former 
is zero since there are as many ways to distribute the 
phases as in (a), as there arc for (b). However, there 
is additional translational entropy for arrangement (a) 
since the domains are twice as large (we thus expect 
AS* = ln2 « 0.69). 

To simulate at fixed m we use Kawasaki dynamics: two 
spins of opposite sign are randomly selected and flipped, 
and the resulting spin conflguration is accepted with the 
Metropolis criterion [2l|. To facilitate frequent transi- 
tions between arrangements (a) and (b) of Fig. [71 we 
also use a collective Monte Carlo move. To this end, wc 
introduce the block domain Bi, which contains all spins 
whose z-coordinate is between Xbi/2 < z < \{bi/2 + 1), 
with bi an integer (periodic boundary conditions must be 
applied). In the collective move, two block domains, Bi 
and B2-, are randomly selected with the constraint that 
|6i — &2I > and even. The domains arc then swapped, 
and the resulting spin configuration is accepted with the 
Metropolis criterion (in our simulations, Kawasaki and 
collective moves are attempted in a ratio 1 : 0.03, respec- 
tively) . 

To test our approach we consider 0.264 < J < 0.278, 
which is between the triple and critical points. We use 
m = —0.426 for this is the value where phases I and III 
were seen to occupy equal volumes. In Fig. |5l^a), we plot 



(a) 


250 




200 




150 




100 




50 







(b) 


640 




160 




40 




10 



ti t\ t\ n 




FIG. 9: Investigation of the critical behavior using a simula- 
tion box with D = 5A; the simulations are performed at fixed 
magnetization m = —0.426 and J = Jcr- (a) The susceptibil- 
ity profiles x{^) for L — 15,20,25,30 (from bottom to top). 
The key point to note from this figure is that x{^) diverges 
with L only at special values z = Zct- (b) Finite-size scaling 
analysis of the average peak height xl of the susceptibility 
profiles. We plot is xl versus L on double- logarithmic scales. 
The dashed line corresponds to a power-law with exponent 
"f/v = 7/4 of the d = 2 Ising model. 
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FIG. 10: The correlation function C{zcr, Az) in the critical 
regime. We show results for D = 5A, L = 30 (a), and D = 
20A, L = 8 (b). The vertical arrow in (a) marks the amplitude 
A of the anti-correlations, which conforms to Eq. ((8|. 



Ini? versus L for J = 0.275; the data are indeed well de- 
scribed by Eq. ([6|), and by fitting 7 can be estimated. In 
Fig.lSKb), we plot the corresponding estimates of 7 versus 
J. Despite the admittedly rather large statistical uncer- 
tainty, our data confirm that the tension is extremely 
low, and that it decreases as J is lowered; both these ob- 
servations are in qualitative agreement with theoretical 
predictions jl8| . 



E. Correlations in the field direction 



In the finite-size scaling analysis of Section IIII CI we 
varied L keeping D = X fixed. We thus assumed the cor- 
relations in the field direction to be short-ranged: critical 
correlations only develop in the lateral L directions, but 
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not in the direction D along the field, such that the result- 
ing critical behavior is effectively two-dimensional (and 
belonging to the d = 2 Ising universality class). To ver- 
ify this assumption we now consider the critical regime 
using a larger value D = 5X. We perform simulations at 
J = Jcr and fixed m = —0.426 (the latter corresponds 
to the average magnetization at the critical point). To 
simulate at fixed m we use Kawasaki dynamics and col- 
lective moves (as in the previous section). However, for 
the collective moves, the block domain was taken to 
be a single lattice layer, containing those spins whose z- 
coordinate equals Zi (at criticality, this choice yields a 
higher accept rate). A pair of layers is chosen randomly 
and swapped, with the constraint that the sign of Vpei{z) 
in the layers is the same, and accepted with the Metropo- 
lis criterion. 

In Fig. HJa), we plot the susceptibility profile x(^) = 
((m(z)^) — (m(2;))^) for various values of L. The sus- 
ceptibility diverges with L only at selected values z = Zcr, 
which "repeat" with the same period as the field. The 
critical behavior is thus spatially confined to those L x L 
slabs for which the corresponding z-coordinate equals one 
of the Zcr- To determine the universality class we com- 
pare the average peak heights xl of the susceptibility 
profiles to the finite-size scaling prediction xl oc L''"/", 
with 7 the susceptibility critical exponent. This result is 
shown in Fig. Mjo), and the d = 2 Ising value j/iy = 7/4 
is strikingly confirmed. Hence, the observed universality 
class does not depend on the value of D used in the scal- 
ing analysis, which a posteriori provides the justification 
for the approach of Section IIII CI 

Next, we ask whether correlations exist between criti- 
cal slabs. To this end, we introduce the pair correlation 
function 

C(zi, Az) cx 

(m(zi) m(zi + Az)) — (m(zi))(m(zi + Az)), (7) 

measured between the slab at z — zi and z — zi + Az, 
respectively. We choose zi to coincide with one of the 
critical slabs, and we normalize such that C(zi,0) = 1. 
In Fig. [TOl we show the correlation function for a system 
with D = 5A (a), and for D = 20 A (b). We find that the 
slabs at Az = n\ with integer n > are anti- correlated 
from the (critical) slab at n = 0. Moreover, the ampli- 
tude A of the anti-correlations is independent of Az, but 
it decreases with D. In fact, an almost perfect "lever 
rule" is observed 



A = A/(D- A). 



(8) 



That is: if there happens to be an excess magnetization 
in one of the critical slabs, the remaining critical slabs 
respond by assuming a lower magnetization, in a manner 
such that the excess magnetization is shared equally on 
average. In the limit D ^ oo, the amplitude A of the 
correlations becomes zero, consistent with our assump- 
tion that long-ranged correlations in the field direction 
are absent. We also point out that the correlations in 
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FIG. 11: The analogue of Fig. [6ja) but for the case d = 2. 
Note the logarithmic vertical scale. The data are obtained 
using fixed D — X. The key thing to note from this figure 
is that the curves for different L do not intersect at a single 
point, implying the absence of a critical point. This, in turn, 
is consistent with d = 1 Ising universality. 



Fig. [TU] are very different from critical correlations; the 
latter decay as power laws, limA2->-oo Cizcr, Az) oc 1/Az'' 
with critical exponent rj, for which we see no evidence in 
our data. In fact, the anti-type correlations of Fig. [TO] 
are also observed in the non-critical regime of the phase 
diagram (explicit checks were performed for J = 0.27 us- 
ing m = —0.82 and m = —0.04, corresponding to a pure 
phase I and phase HI, respectively). 



IV. RESULTS IN d = 2 DIMENSIONS 

We now consider LIC in d = 2 dimensions. The sim- 
ulations are performed on L x D periodic lattices, with 
the field Vpci-{z) again propagating along edge D of the 
lattice. In what follows, the field wavelength A = 8 with 
strength h = 0.1495. 

A. Phase diagram and scaling analysis 

We first determine whether the LIC critical points oc- 
cur in d ~ 2 dimensions also. Since the critical behavior 
was shown to resemble that of a reduced dimension d—l, 
we now expect the universality class of the d = 1 Ising 
model. As is well known, the latter model does not fea- 
ture a critical point. In Fig. 1111 we repeat the cumulant 
analysis of Fig. El^a). In line with the = 1 Ising model, 
we do not observe an intersection point, confirming the 
absence of a critical point. While for small L the curves 
somewhat intersect, the intersections for larger L system- 
atically shift toward larger values of J. Hence, in c? = 2 
dimensions, there is no LIC critical behavior. 

Next, we investigate the fate of the triple point, us- 
ing the same analysis as in Fig. [5] We collect data for 
fixed H = and D = X, while J and L are varied. In 
Fig. I12f a). we plot the Binder cumulant U4 versus J for 
different system sizes L. Consistent with a triple point. 
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FIG. 12: The analogue ol Fig. [S] but for the case d = 2. The 
main difference is that we now observe a critical point, as 
opposed to a triple point, (a) The Binder cumulant as func- 
tion of J for H = and different system sizes L. The curves 
for different L intersect from which we might conclude that a 
phase transition takes place. However, we display in (b) the 
scaling with L of the free energy F{m) at the cumulant in- 
tersection, with F{m — 0) shifted to zero, where we see that 
the depth of the central minimum AF2 — > as L increases, 
while the depth A_Fi of the outer minima appears to be inde- 
pendent of L. This type of scaling is consistent with a critical 
point [13]. 



we observe a sharp intersection, with the value of the 
cumulant at the intersection very close to U4 = 2/3 of 
a triple-peaked distribution. However, the correspond- 
ing free energy is not consistent with a triple point, see 
Fig. HHb), where F{m) is plotted for three different 
system sizes; note that we plot F{m) with the central 
(m = 0) minimum shifted to zero. While F{m) clearly 
reveals three minima, the central minimum does not sur- 
vive in the thermodynamic limit. This can be seen from 
the corresponding "depth", marked AF2 in the figure, 
which decreases with L. In the limit L — ^ 00, we have 
AF2 — )■ 0, and only the outer minima survive, whose cor- 
responding depths then equal AFi. The observation in 
Fig.imjb) that AFi is independent of system size is char- 
acteristic of a continuous transition |27[. Hence, for LIC 
in c? = 2 dimensions, the triple point is destroyed, and 
replaced by a critical point, in this case at Jcr ~ 0.701 
(as expected, this exceeds Jcr, bulk = ln(l-|- V2)/2 of the 
bulk d = 2 Ising model). The LIC phase diagram in 
d = 2 dimensions is thus radically different from d = 3. 
Instead of a "pitchfork" topology, we now have a sin- 
gle line of first-order phase transitions terminating in a 
critical point (Fig. [^Jb)). 




FIG. 13: Free energy F{m) for d = 2 and J = 0.8 > Jcr; the 
data are obtained using D = 2A, H = 0, and two values of L 
as indicated. The key thing to note from this figure is that the 
free energy barrier A.F'i increases with L, indicating a first- 
order phase transition [53]. Note also the spurious minima 
Ml and M2: these reflect meta-stable coexistence states [3] 
whose role in the thermodynamic limit is negligible. 



We still find that, for J > Jcr, the transition is first- 
order. In Fig. [131 we plot the free energy F{m) for 
J = 0.8, iJ = using system sizes L = 40,80 and 
D ~ 2X. The free energy curves are again shifted such 
that F{m = 0) = 0. While for the smaller system the 
minimum at m = is still visible, it has vanished in 
the larger system. In addition, the barrier AJi now in- 
creases profoundly with L, consistent with a first-order 
transition Note also the pronounced fiat region in 

F{m) around to ~ for the larger system: this indi- 
cates two-phase coexistence with negligible interactions 
between the interfaces [s^. When a simulation is per- 
formed in this regime starting from a random initial spin 
configuration, the system phase separates to form struc- 
tures that are strongly affected by the external poten- 
tial; see Fig. [MT a). However, when the system is fully 
equilibrated at a "later time", snapshots show the sys- 
tem containing two coexisting domains of phases I and 
II; see Fig. IWb). In a box with periodic boundaries, 
the domains arrange themselves as two slabs since this 
minimizes the total interface length. 



B. Rounding effects 

Even though the "zebra" phase (phase HI), does not 
survive in the thermodynamic limit in d = 2, we still 
see remnants of this phase in systems of finite size. If 
one simulates at J < Jcr using an appropriate external 
field iJ, one finds that in finite systems F{m) can still 
be cast into the forms of Fig. [ljd,e). Inspection of simu- 
lation snapshots then reveals a condensation of droplets 
onto stripes oriented perpendicular to the field direction 
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(a) (b) (c) 




FIG. 14: Computer generated snapshots obtained using D — lOA and L — 120; white regions correspond to "spin-up". The 
snapshots in (a) and (b) are obtained for J — 0.8 and m = and show "early time" and "late time" snapshots, respectively 
(the simulation was started with a random spin configuration). Since J > Jcr we observe coarsening of domains (a) until two 
large domains have formed (b). In (c), we show a snapshot for the case when J = 0.5 < Jcr and m — 0.52. In this situation, 
domains do not coarsen with time, but remain finite in size, reminiscent of the d = 1 Ising model. 



fFig. [TiT c)). However, the droplet size remains finite in 
this case, owing to the fact that the d = 1 Ising model at 
finite temperature does not support a finite magnetiza- 
tion. Similar finite-size effects occur in colloid-polymer 
mixtures confined to cylindrical pores, which also belong 
to the universality class of the d = 1 Ising model [H, [s^ . 



V. SUMMARY 

In this work, we considered the phase behavior of the 
Ising model exposed to a static periodic field. In = 3 di- 
mensions, we obtain a phase diagram analogous to laser- 
induced condensation observed in fluids undergoing bulk 
liquid-vapor type transitions. That is, a new phase arises 
(the "zebra" phase) and the critical point of the bulk 
model is replaced by two new critical points, and a triple 
point. The main difference compared with fluids is that, 
due to spin reversal symmetry, the corresponding phase 
diagram for the Ising model features a symmetry line. 
The analysis of the present work complements earlier 
works on laser-induced condensation [13, [13 '^^ thed, (i) 
a detailed study of finite-size effects at the triple point 
was presented, (ii) a simple mean-field theory was used 
to elucidate in a qualitative manner how the c? = 3 phase 
transitions depend on the parameters in the external po- 
tential, (iii) a method was presented to measure the ex- 
tremely low tension of interfaces with the zebra phase, 
and (iv) the nature of correlations along the field direc- 
tion was further clarified. 

We additionally considered the fate of laser-induced 
condensation in = 2 dimensions. In this case, we find 
that the zebra phase does not survive in the thermody- 
namic limit, and the corresponding phase diagram fea- 
tures just a single critical point. This critical point oc- 
curs at a temperature below the critical temperature of 
the pure d = 2 Ising model. The universality class of the 
critical point still needs to be determined. The analysis 
of the free energy in Fig. [T^ b) only indicates a critical 



transition, since the barrier AFi is L- independent, but 
no information regarding critical exponents could be ob- 
tained. The practical problem here is that, in computer 
simulations, we are still restricted to system sizes that 
span only a few field wavelengths. We should also men- 
tion that the mean-field theory used in Sec. Ill B predicts 
very similar results in d = 2 as it does in d = 3 and is 
therefore not reliable when applied in d = 2. 

Even though our results were obtained for the rela- 
tively simple Ising model, the generic features of the ob- 
served phase behavior should also apply to real fluids. 
In particular the experimental realization in c? = 2 di- 
mensions should be feasible using a stripe-patterned sub- 
strate. At moderate temperatures, the condensation of 
finite-sized droplets should be observable, while at low 
temperatures a macroscopic demixing should occur (c.f. 
Fig.IUD. 
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